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Context. Sites of massive star formation have complex internal structures. Local heating by young stars and kinematic processes, such as 
Q . outflows and stellar winds, generate large temperature and velocity gradients. Complex cloud structures lead to intricate emission line shapes. 
'*— j CO lines from high mass star forming regions are rarely Gaussian and show often multiple peaks. Furthermore, the line shapes vary significantly 
CZ3 , with the quantum number J up , due to the different probed physical conditions and opacities. 

Aims. The goal of this paper is to show that the complex line shapes of 12 CO and 13 CO in NGC 2024 showing multiple emission and absorption 
features, which vary with rotational quantum number J can be explained consistently with a model, whose temperature and velocity structure 
, are based on the well-established scenario of a PDR and the "Blister model". 
Methods. We present velocity-resolved spectra of seven 12 CO and 13 CO lines ranging from J up = 3 to J up = 13. We combined these data with 
f**«« | l2 CO high-frequency data from the ISO satellite and analyzed the full set of CO lines using an escape probability code and a one-dimensional 

■ full radiative transfer code. 

Results. We find that the bulk of the molecular cloud associated with NGC 2024 consists of warm (75 K) and dense (9 • 10 5 cm 3 ) gas. An 
additional hot (~ 300 K) component, located at the interface of the HII region and the molecular cloud, is needed to explain the emission of the 

i—{ • high-J CO lines. Deep absorption notches indicate that very cold material (~ 20 K) exists in front of the warm material, too. 

Conclusions. A temperature and column density structure consistent with those predicted by PDR models, combined with the velocity structure 

0^ . of a "Blister model", appropriately describes the observed emission line profiles of this massive star forming region. This case study of NGC 

■ 2024 shows that, with physical insights into these complex regions and careful modeling, multi-line observations of I2 CO and 13 CO can be 
' used to derive detailed physical conditions in massive star forming regions. 
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1 . Introduction The massive star forming region we selected for our study 

is the HII region NGC 2024 and its associated molecular 

High-mass star forming regions are very complex. The inter- , , u - . ■ i . j » j- . t , 1C , » 

b b b cloud, which is located at a distance of 415 pc (Anthony- 
action of OB stars, which are often deeply embedded, with m i lrioo i T „ . . £ „ . n „ . , . . 

r J Twarog 1982). It is a part of Orion B, a well-studied giant 

the surrounding molecular cloud via radiation and outflows molecular doud (e Maddalena et aL[M Lada et al. [19911 

and condensations of cold gas, still largely unaffected by the ^ ^ rjggg Mhchell et a,|20m 
star forming activities, lead to a complex density and tern- 

perature structure, causing intricate shapes of the observed A P ossible scenario for the three-dimensional structure of 

emission lines. Some CO lines in such regions (e.g., M 17, the NGC 2024 re 8 lon was Prosed by Barnes et al. ([1983, 

Stutzki & Guesten et al. M W3, Kramer et al. [2004] and who combined observations from optical to radio wavelengths 

Mon R2, Giannakopoulou et al. \mj} have multiple peaks, to construct a model of the cloud ( Fig. 0. In their model, the 

most likely due to absorption by foreground material. The line HI1 re g lon slts ln front of the bulk of the cloud but 18 P 3 ^ ob " 

shapes vary significantly, depending on the rotational level and scured ln the °P tlcal b y a ver y Prominent ridge of cold molec- 

the observed isotopes, which trace different regimes of physical ular material. This geometry provides an explanation for the 

conditions and optical depths. complex line shapes, which may consist of contributions from 

fore- and background regions, including self-absorption. The 

Send offprint requests to: M. Emprechtinger ionizing sources of NGC 2024 are invisible at optical wave- 

* Current address: California Institute of Technology, Owens Valley lengths because they are obscured by the dust ridge. IRS2 b, a 

Radio Observatory, Big Pine, CA 93513 USA late O or early B star, has been identified as the primary ioniz- 
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Fig. 1. Schematic view of the geometry of NGC2024. 
Abbreviations are as follows: DMC = dense molecular cloud, 
EMC = extended molecular cloud, lF=ionization front and t 
= high optical depth cloud (optical dust bar). Adapted from 
Barnes et al. ( [T989l > by Giannini et al. (I20001 . 



ing source of NGC 2024 through observations in the near in- 
frared (Bik et al. 2003). The extinction (A v ) through the dust 
ridge along the lines-of-sight to stars inside the HII region are 
in the range of 15 to 25 mflS . Furthermore, the HII region is ex- 
pected to expand into the molecular cloud and to trigger star 
formation. Indeed some protostellar condensations (FIR 1 to 7) 
have been detected close to the HII region (Mezger et al. 1988 
& 1992). Their masses, derived from sub-millimeter continuum 
emission, range from 1.6 to 5.1 M (Visser et al 1998). Because 
these masses correspond to visual extinctions between 21Q mag 
and 870"" w and the optical depth of the dust ridge at the near 
side of the H II region is only 15-25'" flg , it is very likely that 
these condensations are located in the dense molecular cloud 
(DMC) behind the HII region. Massive outflows have been de- 
tected close to the source FIR 5 (e.g., Richer 1990, Sanders & 
Willner 1985 ) giving additional evidence for ongoing star for- 
mation. 

Detailed investigations of the molecular cloud associated 
with NGC 2024 were performed by Graf et al. ([19901 [19931 . 
who studied multiple 12 CO lines (up to J=7-6) and its isotopo- 
logues. They derived temperatures of 23.5 K and 67.4 K for 
the foreground and the background component, as defined by 
Barnes et al. (1989), respectively. The corresponding column 
densities are 5.4 • 10 22 ctrT 2 (A v = 56 mag ) and 2.0 ■ 10 23 cirT 2 
(A v = 210 mag )). Furthermore they found different velocities 
of 9 km/s and 1 1 km/s for the foreground and the background 
component, respectively. The peak velocities of the main com- 
ponent of the 12 CO lines, Vlsr =13 km/s, differ significantly 
from those of the optically thin lines (e.g., of C ls O), which are 
at a velocity of 11 km/s. Graf et al. (1993 ) suggested that the 
blue side of the 12 CO lines are absorbed by the dust ridge at a 
velocity of 9 km/s. 



This velocity shift between background and foreground 
component can be explained by a "Blister Model" (Israel[l978 
Zuckerman 1 1 9731 l. In such a model, an OB star is assumed to 
be located close to the surface of the molecular cloud form- 
ing an HII region. The ionization front of such an HII region 
moves slowly into this cloud establishing a high pressure gra- 
dient. Because of this pressure gradient the ionized gas moves 
away from the molecular cloud. Therefore hydrogen recombi- 
nation lines appear at negative line of sight velocity offsets of 
about 3 km/s relative to the molecular lines, if the HII region 
is located at the near side of the molecular cloud (Israel 1978). 
In the case of NGC 2024 the ionized material cannot flow into 
space, but pushes on the foreground material. We therefore ex- 
pect that the foreground component is at somewhat lower ve- 
locities than the bulk of the molecular material. The assumption 
that the velocity structure of NGC 2024 is indeed caused by a 
"Blister" is buttressed by the fact that the H109o' recombina- 
tion line appears at 7 km/s (Israel [19781 1. and thus even more 
blueshifted. 

C is O 2-1 and C 17 2-1 emission observed by Graf et 
al. (Q993 ), follows the 1 .3 mm continuum very well and peaks 
close to the far infrared sources detected by Mezger et al. ( 1988 
& [19921 . Contrary to the C ls O and C n O lines, the map of the 
integrated intensity of optically thicker 12 CO 7-6 line shows 
features similar to those seen in the 6 cm continuum map 
(Crutcher et al. 1986), which traces the ionized gas. This in- 
dicates that the 12 CO 7-6 mainly originates from a photo dom- 
inated region (PDR) at the surface of the molecular cloud and 
not from embedded protostellar objects. The excitation condi- 
tions in this PDR were studied using the integrated intensity 
of far infrared line emission observed with the ISO satellite 
(Giannini et al 120001 . They studied [Nil], [NIII], and [OIII] 
lines from the HII region and high- J CO lines, [CII], and [OI], 
which originate (in part) from the PDR. Comparison with the 
predicted line ratios of a PDR model (Burton et al. 199(3 re- 
vealed a density of « 10 6 cm 3 and a UV field of 3 • 10 4 G at 
the surface of the molecular cloud. 

PDR models for high densities and high UV fields, as rel- 
evant for NGC 2024, predict a column density for hot CO 
(> 100 K) of the order of 10 16 - 10 17 cirT 2 (see, Rollig et 



al.|200f]). High-J CO lines (in this paper we refer to all lines 
with frequencies > 1 THz as high- J) would be emitted exclu- 
sively from such hot and dense material, because of their high 
critical density (~ 10 7 cirT 2 ) and high energy of the upper level 
(E up > 250K). 

In this paper we combine twelve transitions of 12 CO and 
13 CO, ranging from J up =3 to J up = 19 to trace molecular mate- 
rial NGC 2024 under very different physical conditions. The 
goal of this paper is to find, based on the above mentioned 
scenarios ("Blister Model" and PDR) a model for the veloc- 
ity, density and temperature structure of NGC 2024, which re- 
produces the observed line intensities and line shapes consis- 
tently. For this purpose we observed seven CO lines spectrally 
resolved which are presented in section[3] In sectional we first 
examine the physical conditions of the gas, which emits the 



1 All of their results are available 
http://www.phl.uni-koeln.de/pdr-comparison 
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high-J CO lines, using an escape probability code and compare 
the result with the expectation from PDR models. The outcome 
of the escape probability calculations are used as input of a five 
component radiative transfer model (section l4~2b . An interpre- 
tation of the model is given in section [5] 

2. Observations 

2.1. Observations of 12 CO J = 13 - 12 at APEX 

The 12 CO 13-12 (v = 1496.922909 GHz, Miiller et al. EOOfl ) 
observations were carried out on November 22, 2005 at 
the APEX 12 m telescope (Giisten et al. 2006), located on 
Llano de Chajnantor, Chile. We used the CO N + Deuterium 
Observations Receiver (CONDOR, Wiedner et al. |20061 l, a 
heterodyne receiver that operates at THz frequencies (1.25- 
1.53 THz). The typical double sideband (DSB) receiver noise 
temperature was between 1500 K and 1900 K. The mean 
atmospheric transmission at zenith during the observations 
was ~ 0.2. As a backend, we used the APEX Fast Fourier 
Transform Spectrometer (FFTS), which has an intrinsic resolu- 
tion of 60 kHz. 

The expected diffraction limited main beam size of a 12 m 
telescope at 1.5 THz is 4.3". Because the accuracy of the pri- 
mary surface derived from planet observations at frequencies 
between 350 GHz a nd 1500 GHz (with CONDOR) is 18 fim 
(Giisten et al. 2006), we expect some of the power to be di- 
rected into an error beam. We assume that the error beam is ap- 
proximately 80", which is derived from the approximate size of 
the individual panels of the telescope (~ 70 cm). To determine 
the efficiency for the calibration to the main beam tempera- 
ture (Tmb), we measured the continuum of Moon and Mars 
and compared the observed values with models. These mea- 
surements lead to beam efficiencies of 0.4 and 0.09-0. 1 1 for 
the Moon and Mars, respectively. The different beam efficien- 
cies arise because Mars is about one-fourth the size of the error 
beam (18.2" on the date of the observations), and thus these 
observations suffer from beam dilution. A more detailed dis- 
cussion of the beam sizes and coupling efficiencies is given in 
Volgenau et al. (120081) and Wiedner et al. J2006I) . 

Since NGC 2024 is an extended source, clearly larger 
than the error beam (Graf et al. 119931 Kramer et al. 119961 ), 
we use a beam efficiency of 0.4 for the calibration to Tmb- 
Pointing and focusing were determined from observations of 
Mars, and we estimate a pointing accuracy better than 7". We 
observed 12 CO J=13-12 towards five positions, along a line 
from FIR 5 via IRS 3 to IRS 2. The coordinates of FIR 5 
are a = 5 h 41 m 44.18 s , 6 = -1° 55' 38.0" (J=2000, Mezger 
et al. 1 1988b . The offsets of the other positions with respect to 
FIR 5 are listed in Tab. □ 

2.2. Observations of mid-J CO lines with KOSMA 

In addition to the 12 CO 13-12 transitions, we used archival, 
so far unpublished, maps of I2 CO 6-5 (v = 691.5 GHz), 
12 CO 3-2 (v = 345.7 GHz), 13 CO 6-5 (v = 661.1 GHz) 
and 13 CO 3-2 (v = 330.6 GHz). The observations were con- 
ducted between January 27 and February 2 1998 at the Kolner 



Table 1. Positions of the 12 CO J=13-12 observations. The (0,0) 
position is a = 5 h 41 m 44.18 s , 6 = -1° 55' 38.0" (J=2000, 
Mezger et al. [19881 . 



Position 


Aa ["] 


A<5 ["] 


FIR 5 


0.0 


0.0 


IRS 3 


3.3 


16.0 


#1 


8.8 


28.6 


#2 


13.8 


41.9 


IRS 2 


24.6 


68.5 



Observatorium fur Submillimeter Astronomie (KOSMA), lo- 
cated on Gornergrat, Switzerland. The beam sizes (HPBW) of 
these observations are ~ 50" and 82" at 690 GHz and 345 GHz, 
respectively. To map the central part of the molecular cloud as- 
sociated with NGC 2024, we used a dual-channel SIS receiver 
(Graf et al. [1998 ) with typical DSB receiver noise tempera- 
tures of 100 and 400 K for the J=3-2 and J=6-5 lines, respec- 
tively. As backends, we used two acousto optical spectrome- 
ters (AOS). The Medium Resolution Spectrometer (MRS) has 
a bandwidth of 0.3 GHz and a resolution of 360 kHz. The Low 
Resolution Spectrometer (LRS) has a bandwidth of 1 GHz and 
a resolution of 1150 kHz. The forward efficiency, F e g = 0.93, 
was determined by skydips. For the J=3-2 observations, we 
used beam efficiencies (B e ff) of 0.59 and 0.62 for the 12 CO 
and the 13 CO line, respectively. For the J=6-5 observations, the 
beam efficiencies were 0.40 and 0.48, respectively. 

2.3. Observations of mid-J CO lines with NANTEN2 

Simultaneous observations of 12 CO 7-6 (v = 806.65 GHz) 
and 12 CO 4-3 (v = 461.04 GHz) were carried out with the 
NANTEN2 4m telescope at Pampa La Bola, Chile. We ob- 
tained a 2' x 2' map, which was centered on NGC 2024 IRS 3. 
The observations were carried out in December 2007, using 
the dual-channel 460/810 GHz test receiver, which had DSB 
receiver temperatures of ~ 750 K and ~ 250 K for the up- 
per and the lower channel, respectively. Two AOS with band- 
widths of 1 GHz and channel widths of ~ 560 kHz were used 
as backends. The beam sizes (HPBW) of the observations were 
25" and 37" for the 12 CO 7-6 and 12 CO 4-3 observations, re- 
spectively. The beam efficiencies are 0.5 and 0.45 for 460 and 
810 GHz, respectively, and a forward efficiency of 0.86 was 
measured at both frequencies. The position of IRS 3 was taken 
to be a = 5 h 41 m 44.40 s , 6 = -1° 55' 22.8" (J=2000, Mezger 
et al. 119921 ). Pointing was checked on IRc2 in Orion A right 
before the observations and is expected to have an accuracy of 
< 7". 

3. Observational Results 

To be able to compare the observed lines with each other, we 
convolved the maps observed with the KOSMA telescope and 
the NANTEN2 telescope to a spatial resolution of 80", which is 
the spatial resolution of the J=3-2 spectra as well as the approx- 
imate size of the error beam of the 12 CO 13-12 observations. 
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Since 70%-80% of the radiation is expected to come from the 
error beam at these high frequencies, we considered 80" as the 
spatial resolution of the 12 CO 13-12 observations as well. An 
analysis of the maps of the lower-J lines will be given in a sub- 
sequent paper. 

The spectra are shown in Fig. [2] and the line parameters are 
listed in Tab. [2] While most of the lines have a Gaussian line 
profile ( 12 CO 13-12, 13 CO 6-5) or only relatively weak blue 
shifted shoulders ( 12 CO 7-6, 12 CO 6-5, 13 CO 3-2), the 12 CO 
4-3 and 3-2 lines show complex line shapes with absorption 
notches, enhanced emission from the red shifted wing, and an 
additional bump at ~ 4.5 km/s. 

The emission of the J=3-2 and J=6-5 transitions of 12 CO 
and 13 CO are fairly uniform in line intensity and shape at all 
five positions, and the integrated intensity of these lines drops 
towards IRS 2 only by ~ 30%. This uniformity is expected, 
since the separation between the positions of FIR 5 and IRS 2 
is only 73", which is on the order of the spatial resolution of the 
observations. The 12 CO 13-12 line peaks towards the southern 
three positions (FIR 5, IRS 3 and position #1) and declines no- 
ticeably towards the north, indicating that high-J CO emission 
originates from a rather compact region around IRS 3. 

Because of the different and complex line shapes of the 
12 CO and 13 CO lines, we give the centroid velocity of the 
spectra in the following. The velocities of the 13 CO lines are 
~ 11 km/s, whereas the velocity of 12 CO 3-2 and 12 CO 6-5 
is ~ 12 km/s. The J=13-12 line of 12 CO is found at a veloc- 
ity of ~ 13 km/s, although there is a variation in the veloc- 
ity of about 1 km/s with position. The velocity of the 13 CO 
lines is uniform throughout all five positions, whereas the ve- 
locity of 12 CO 3-2 and 12 CO 6-5 declines from south to north. 
This is due to a redshifted outflow detected in these two lines, 
which is strongest at the southern positions and reaches veloc- 
ities up to ~ 20 km/s. This outflow was detected previously 
by Sanders & Willner (1985). The lack of a blueshifted wing 
in our spectra is consistent with previous observations (e.g., 
Richer 1990 & Richer et al. 1992). One explanation given by 
Richer et al. (1989) is that the outflow exists in a very complex 
region close to the ionization front of the HII region. North of 
its origin lies very dense gas seen in 1.3 mm map (Mezger et 
al. [19921 1 and HCO + (Richer et al. |T989b . Any material ejected 
by the protostar in this direction might be retarded by dense gas 
and/or destroyed by strong UV-radiation. 

The 12 CO 13-12 data displayed in Fig. [2] are binned to a 
resolution of 0.48 km/s, which is sufficiently high, because the 
line width of the J=13-12 line is on the order of 2.5 km/s. At 
this resolution we detected the line in four out of five positions 
with an S/N ratio > 3<x. At the position of IRS 2, we can claim 
only a tentative detection of 7.5 + 3 K. 

The 12 CO 6-5 line is not symmetric, but shows a shoul- 
der at its lower velocity side. Its velocity of ~ 12 km/s is ap- 
proximately 1 km/s redshifted with respect to the 13 CO 6-5 
line. The difference in velocity can be explained by the two- 
component model, suggested by previous studies (Barnes et 
al.[l989, Graf et al. 1 1993b . In this model, the foreground com- 
ponent (r in Fig. Q]) is assumed to be at a lower Vlsr than 
the background component (9.2 km/s and 11.1 km/s, respec- 
tively). Furthermore, the temperature of the foreground corn- 



Table 2. Here the integrated intensities, Vlsr and Av of the lines 
shown in Fig.|2]are listed. Because of the complex line shapes, 
especially from the lower-J 12 CO lines, the Vlsr and Av have 
been determined by using the first and second moment of the 
spectra, respectively. 



Position 


1 T ra b 


dv 


Vlsr 


Av 




[Kkm/s] 


[km/s] 


[km/s] 


12 CO 13 - 12 


FTP ^ 


120 ± 


10 


1Z.ZJ + U.J 


J.J^f ± U.J 


IRS 3 


109 ± 


10 


12.75 ± 0.3 


3.92 ± 0.5 


# 1 


127 ± 10 


13.50 + 0.3 


4.51+0.5 


#2 


57 ± 


10 


13.17 + 0.3 


3.54 ± 0.5 


IRS 2 


< 30 




_ 






U. 


C07-6 




TRS 


229 ± 


10 


11.5 + 0.2 


4 ^ + S 






12, 


C06-5 




FTP S 
FlJX J 


Too i 

Zoo ± 


10 


12.0 + 0.15 


7 s -t- n s 

/ . J ± U.J 


IRS 3 


289 + 


10 


11.9 + 0.15 


7.2 ± 0.5 


# 1 


279 + 


10 


11.9 + 0.15 


7.0 + 0.5 


# 2 


265 ± 


10 


11.9 + 0.15 


7.2 + 0.5 


IRS 2 


235 ± 


10 


11.8 + 0.15 


7.5 + 0.5 








C06-5 




rlK j 


159 i 


:5 


10.9 + 0.15 


Z.S ± U.J 


IRS 3 


157i 


:5 


10.9 + 0.15 


2.7 + 0.5 


# 1 


149 i 


:5 


10.8 + 0.15 


2.6 + 0.5 


#2 


138i 


:5 


10.7 + 0.15 


2.5+0.5 


IRS 2 


106 4 


:5 


10.5 + 0.15 


2.3 + 0.5 






12, 


CQ4-3 




IRS 3 


257.8 i 


t 10 11.5 + 0.2 


11.8 + 0.3 






12, 


C03-2 




FIR 5 


213 ± 


10 


11.9 + 0.15 


10.2 + 0.2 


IRS 3 


202 ± 


10 


11.9 + 0.15 


11.3 + 0.2 


# 1 


186 + 


10 


11.6 + 0.15 


10.7 ± 0.2 


#2 


171 ± 


10 


11.0 + 0.15 


8.55 + 0.2 


IRS 2 


155 ± 


10 


10.7 + 0.15 


6.52 ± 0.2 






13, 


CO 3 -2 




FIR 5 


146 4 


:5 


10.9 + 0.1 


3.36 + 0.2 


IRS 3 


144 4 


:5 


10.8 + 0.1 


3.42 ± 0.2 


# 1 


136 4 


:5 


10.8 + 0.1 


3.34 + 0.2 


#2 


127 4 


:5 


10.8 + 0.1 


3.24 + 0.2 


IRS 2 


109 4 


:5 


10.7 + 0.1 


2.93 + 0.2 



ponent (24 K) is lower than the temperature of 67 K of the 
background (Graf et al. 1993). Thus, the foreground component 
adds an absorption feature centered at 9.2 km/s to the 12 CO 6- 
5 line emitted from the background component. Therefore, the 
centroid velocity of the spectra appears redshifted with respect 
to the background component. In 13 CO 6-5, no signs of self- 
absorption appear, and thus its centroid velocity coincides with 
the Vlsr of the background component. 

The profile of the 13 CO 3-2 line looks similar to the 
12 CO 6-5 line, but it does not appear redshifted and lies at a 
velocity of ~ 1 1 km/s. The shape of the 12 CO 3-2 lines is dom- 
inated by deep absorption notches at velocities of 9.8 km/s and 
11.9 km/s, which are caused by material located in front of the 
HII region. The alternative scenario that the emission is com- 
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Fig. 2. Spectra of CO emission from pointed observations in NGC 2024. The spectra are all convolved to a resolution of 80", 
which is approximately the resolution of the observed 12 CO 3-2 and the 13 CO 3-2 spectra as well as the resolution of the error 
beam of the 12 CO 13-12 observations. The spectra are taken towards the positions (north to south left to right) as denoted in 
Tab.Q] The grey spectra, overlaid on I2 CO 6-5 and 12 CO 3-2 at ISR 3, are the 12 CO 7-6 and the 12 CO 4-3 spectra, respectively. 



posed of three individual, equally strong (~ 25 K), velocity 
shifted cloud components is unlikely, because these compo- 
nents are not seen in most of the other lines. A fourth, weaker 
peak (~ 3 K) can be seen at 4.6 km/s. I2 CO 2-1 observations 
(Kramer et al. 1996) revealed that this component extends fur- 
ther to the north-east and seems to be kinematically distinct 
from the main component. 

The spectra of 12 CO 7-6 and 12 CO 4-3, convolved to a res- 
olution of 80", are superimposed on the 12 CO 6-5 and 12 CO 3- 
2 in Fig. [2] respectively. The 12 CO 4-3 spectrum looks very 
much like the one of 12 CO 3-2. The two absorption notches, 
the 4.6 km/s component, and the redshifted outflow can be seen 
in both spectra. Even the intensities are nearly identical. The 
only difference is that the emission from the outflow, i.e., at 
Vlsr > 12 km/s, is slightly stronger in 12 CO 4-3. 

The blueshifted shoulder of the 12 CO 6-5 spectrum, inter- 
preted as self absorption, appears similar but more pronounced 
in 12 CO 7-6. However, the 12 CO 7-6 line is narrower, because 
the outflow signal is weaker than in 12 CO 6-5. 12 CO 7-6 seems 
to suffer much more from absorption and what is seen as a 
shoulder in 12 CO 6-5 appears as a second peak here. The dif- 
ferent strengths of the absorption of I2 CO 7-6 and 12 CO 6-5 
requires a strong gradient in the excitation of those lines. The 



velocity of the 12 CO 7-6 (1 1.5 km/s) line is in between the ve- 
locity of 12 CO 6-5 (~ 12 km/s) and 13 CO 6-5 (~ 11 km/s). 
If these velocity differences are indeed due to foreground ab- 
sorption, the velocity dispersion of the 12 CO 7-6 absorbing ma- 
terial is significantly lower than the velocity dispersion of the 
material which absorbs 12 CO 6-5. The minimum, most likely 
caused by the foreground absorption, lies at 9.5 km/s, which 
is 0.3 km/s blue-shifted with respect to the absorption notches 
seen in 12 CO 4-3 and 12 CO 3-2. 



Graf et al. ([T9931 observed NGC 2024 in 12 CO 7-6 as well, 
using the UKIRT telescope, which has a similar spatial reso- 
lution. Their spectrum at FIR 5, 16" south of IRS 3, shows a 
Vlsr of 12.8 km/s, clearly red shifted with respect to our obser- 
vations. However the frequency stability of the Laser system 
that was used as local oscillator in these early measurements 
was about 1 MHz (~ 0.4 km/s). Furthermore, the UKIRT ob- 
servations were made in a double beam switch mode with a 
chop throw of 3'. Thus their off-position might have been not 
completely clean, leading to a different apparent line shape. 
Possible pointing errors cannot be the reason for different ve- 
locities reported in Graf et al. ( 119931 1 and this work, because no 
big velocity gradients with position are seen in NGC 2024. 
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4. Modeling Results 

The aim of this investigation is to show that the major fea- 
ture of the complex shapes and intensities of all observed I2 CO 
and 13 CO lines can be explained consistently with a physically 
plausible scenario, which is based on the "Blister model" and 
the PDR scenario (see section[TJi. 

4.1. Escape Probability Code Results 

In this section we examine the physical properties of the 
high-J CO lines using an escape probability code (Stutzki & 
Winnewisser 1985 ). By comparing the column density of hot 
CO required to explain the observed high-J CO emission with 
the column density expected by PDR models we check, if the 
assumption of a PDR-scenario is consistent. Furthermore we 
use the results of this models as a first guess for a more sophis- 
ticated five component model (section POK which reduces the 
free parameter of this model. 

In addition to the new I2 CO 13-12 spectra of NGC 2024 
presented in this work, the integrated CO intensities from 
Jup= 14 to J up = 17 observed with the ISO satellite towards FIR 5 
(Giannini et al. 2000) were taken from literature. The 12 CO 13- 
12 line does not show any sign of absorption or emission at 
~ 9, 5 km/s. Thus, we assume that the emission of I2 CO 13-12, 
as well as of the other high-J CO lines, is purely determined by 
the material located behind the HII region. Therefore, we can 
use an escape probability code (Stutzki & Winnewisser 1985 ) 
to model the emission of 12 CO J > 13. In this code, the emit- 
ting gas is treated as an isothermal cloud with homogeneous 
density. For a given set of kinetic temperatures (Tkj n ), H2 den- 
sities (n(H_2)), and 12 CO column densities per velocity interval 
(N(CO)/Av), the emitted line strengths (main beam tempera- 
tures and integrated intensities), as well as the optical depths 
(t) at the line center are calculated. 

To match the observed intensities of the high-J lines, a 12 CO 
column density of 9.5 + 0.7 • 10 16 crrT 2 in a velocity interval 
Av = 4 km/s, i.e., the width of the 12 CO 13-12 line, is re- 
quired. This column density is in good agreement with PDR 
models (see section[TJ. Constraints on the gas temperature and 
H.2-density, however, are looser; many solutions are possible. 
The two best fit solutions (both with a x 2 of 2.9) yield a H2- 
density of 1.3 ■ 10 6 cirT 3 and a temperature of 250 K and 
nn 2 = 3 ■ 10 5 cm 4 and T=410 K, respectively. All models with 
ax 2 lower than 10 have in common that T 3 • n is approximately 
constant. 

If we assume that the gas is distributed uniformly through- 
out the (80") beam (and the area filling factor is one), then this 
layer of the PDR has a thickness of 5 • 10 14 to 1.5 • 10 15 cm (= 
33-100 AU). At a distance of 415 pc, 100 AU corresponds to 
0.25", which indicates that this hot component is indeed a thin 
layer at the, possibly clumpy, surface of the molecular cloud. 

4.2. Full Radiative Transfer Model 

To get a more complete picture of the source, which explains 
also the complex line shapes of the lower-J lines we have to 
use a multilayer radiative transfer model, which treats multiple 



emission components and self-absorption correctly. With such 
a radiative transfer model, we fitted the line shapes and inten- 
sities of all seven 12 CO and 13 CO lines observed at the position 
of IRS 3 and the integrated intensities of 12 CO with J up > 13, 
simultaneously. The I2 CO line with J up > 13 are the same data 
we used in the previous section. These data have been observed 
at the position of FIR 5, which lies 16" south of IRS 3, a separa- 
tion much lower than the spatial resolution of our observations. 
We used SimLine (Ossenkopf et al. 2001), a 1 -dimensional, 
spherical radiative transfer code. SimLine computes the pro- 
files of molecular rotational lines for an arbitrary density, tem- 
perature and velocity distribution, specified as a set of discrete 
layers, by integrating the radiative transfer equation numeri- 
cally. The population of the individual levels of a molecule are 
computed by solving the full system of balance equations itera- 
tively. By setting the inner radius of the model to a value, which 
is much larger than the beamsize times the distance of the ob- 
servations, we used SimLine in a quasi plane-parallel way. 

4.2.1 . Constraints on the models 

A close look at the data in Fig. [2] and Table [2] reveals an ap- 
parent contradiction in the line strengths at v L sr =11 km/s. 
The strong (~ 50 K) 13 CO 6-5 emission, which indicates a 
large column density of hot molecular gas, belies the relatively 
weak (< 15 K) 12 CO 13-12 emission at the same velocity. 
The possibility that the 12 CO 13-12 line is absorbed by the 
foreground gas can be ruled out, since the estimated column 
density (Nh, = 2 ■ 10 22 cirT 2 , which corresponds to A v =20- 
25" Mg , Bik et al. [2003 ) is orders of magnitudes too low for 
the derived temperatures (~ 25 K) to cause this absorption 
(Graf et al. 119931 ). Relatively strong 12 CO 13-12 emission at 
a vlsr of ~ 13 km/s, where no other optically thin line shows 
an emission peak, suggests an error in the frequency calibra- 
tion. The origin of such an error is puzzling, given that par- 
ticular care was taken to check calibration during the observa- 
tions. Immediately prior to observing NGC 2024, we observed 
Orion FIR 4 to check telescope pointing, and the 12 CO 13- 
12 emission had the expected velocity (Wiedner et al. 120061 
Kawamura 120021 Wilson et al. [2001). The uniformity of the 
Vlsr throughout the 13 CO 6-5 map also diminishes the likeli- 
hood that the anomalous 12 CO 13-12 velocity is due to a point- 
ing error. Moreover, all software-set values (e.g., rest frequency 
and sky frequency of the line) were checked carefully, but no 
error was found. Despite this scrutiny we are forced to assume 
that the v LSR of 12 CO 13-12 and 13 CO 6-5 are both 10.9 km/s, 
to get a reasonable solution for our model. This assumption is 
buttressed by unpublished high-J CO observations (Marrone, 
priv. comm.), which also indicate a Vlsr of ~ 1 1 km/s. 

In addition to the previous two component models we in- 
troduced a hot, dense, and thin gas component, located at the 
interface of the HII region and the DMC (Fig. [TJ. An indepen- 
dent estimation of the physical conditions of this hot material, 
which is illuminated by OB stars within the HII region, has 
been derived with the escape probability code. In the case of the 
foreground material, we have to assume several components as 
well. First, we have to adopt two components of cold gas at dif- 
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Table 3. Physical properties of a 5-component radiative transfer model of NGC 2024. Properties with a range of values vary 
continuously with depth. 



Component 


T[K] 


n(H 2 ) [cm" 3 ] 


r [cm] 


N(H 2 ) [cm" 2 ] 


VLSR[km/s] 


Avnu-b [km/s] 


Bl 


75 


9- 10 5 


8.0- 10 16 


7.2 ■ 10 22 


11.0 


1.8 


B2 


75-330 


9 • 10 5 -2 • 10 6 


6.2- 10 13 


8.9- 10" 


11.0 


1.8 


Fl 


330-40 


1 • 10 5 


8.9 • 10 14 


8.9 • 10" 


9.3 


1.3 


F2 


40-30 


1 • 10 5 - 3 • 10 4 


1.7 • 10 17 


1.04- 10 22 


9.3-9.7 


1.3-1.1 


F3 


30-15 


5 • 10 4 


5.9- 10 15 


2.9 • 10 20 


12.1 


1.8 



ferent velocities, as both the I2 CO 4-3 and the 12 CO 3-2 spectra 
show two absorption notches. In addition, it is reasonable that 
material located at the interface of the foreground component 
and the ionized gas is heated as well, and thus we introduced a 
layer of hot gas on this side of the HII region, too. Furthermore, 
the total column density of all components located in front of 
the HII region is restricted by the measurements of Av towards 
the obscured stars, which translates into a total H2 column den- 
sity of 2 • 10 22 cm" 2 . 

4.2.2. Best Fit Model: Reproducing the observed line 
profiles 

The result of the modeling is listed in Tab. [3] As mentioned 
above, five components (three foreground and two background 
components) were appropriate to get a reasonable solution for 
the observed spectra. Each these five components is described 
by 5 parameters (7^ nffy), r, vi$r and Av) and thus the model 
contains 25 parameters. However, since the hot foreground 
component (Fl) is determined by the temperature and column 
density of the hot background component (B2) and the density, 
Vlsr and Av of the bulk of the foreground material (F2), these 
five parameters are not free. In addition, we set Vlsr and Av of 
the background components (Bl and B2) to the same values, 
and thus the number of free parameters is 18. But even these 
18 parameter are not fully free. The total column density of 
the foreground components (F1+F2+F3)) must correspond to 
an A v of \Q-25 mag to match the IR-observation, the tempera- 
ture and the column density of the hot background component 
(B2) has to be consistent with the prediction of the PDR mod- 
els, and the velocity difference between the background com- 
ponents and the bulk of the foreground material must be neg- 
ative and in the range of a few km/s in order to be explained 
by the blister model. Therefore the number of effectively free 
parameters is ~ 15. 

The background material was divided into two subcompo- 
nents, denoted by Bl and B2. Bl consists of warm (T= 75 K) 
and dense gas (n(H2)= 9 • 10 5 cm -3 ) with a high column den- 
sity (NXH2) = 7.2 • 10 22 cm" 2 ). B2 is a thin layer with a steep 
temperature and density gradient. At the surface of the cloud 
the temperature reaches 330 K and the density is as high as 
2 • 10 6 cm" 3 . Both Bl and B2 have a velocity of Vlsr = 11 km/s 
and a FWHM of the turbulent velocity of Av tur b = 1.8 km/s. 

The foreground component was divided into three subcom- 
ponents, Fl, F2, and F3, to explain the two absorption dips in 
the 12 CO 4-3 and 12 CO 3-2 spectra and for physical reasons, 



as mentioned above. The Vlsr of F2 shows a gradient from 
9.3 km/s to 9.7 km/s towards the observer, and the Vlsr of F3 
is 12 km/s. The total column density of the foreground com- 
ponent (F1+F2+F3; N(H 2 ) to t = 108 • 10 22 cm" 2 ) correspond 
to an optical extinction of ~ ll mag , which is slightly lower 
than the 15 to 25 mag found by Bik et al. (120031 . To convert 
the column density into A v we used the relation Nh, [cm" 2 ] = 
0.94 x 10 21 A V (mag). However, since the light of some of these 
stars might be extincted by dust in the vicinity of the star, the 
lower value of A v is probably more appropriate. 
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Fig. 3. Integrated line intensity of 12 CO versus J up . The ob- 
served data (full circles) include all 12 CO lines from this 
work and observations with the ISO satellite from Giannini et 
al. (f2000). The solid line shows the results of the full radiative 
transfer code SimLine 

The complex shapes of the emission lines, especially the 
deep absorption notches in 12 CO 3-2, require depth-dependent 
gradients in the foreground subcomponents. F2 contains most 
of the foreground material (> 96 % of the mass). It consists 
of a 1.54 • 10 17 cm deep, 40 K warm gas layer, whose density 
decreases linearly from 10 5 cm" 3 to 3 ■ 10 4 cm" 3 . Subsequently 
the F2 gas density stays constant, but the temperature drops 
to 30 K within the next 1.4 • 10 16 cm. Simultaneously to the 
temperature decrease, the Vlsr shifts from 9.3 km/s to 9.7 km/s. 
F3 is a comparatively thin (Nh 2 = 2.9 • 10 20 cm -2 ) component 
with a constant density of 5 • 10 4 cm" 3 and a temperature that 
drops from 30 K to 15 K. The Vlsr of F3 is 12.2 km/s. F3 is 
responsible for one of the two absorption dips in I2 CO 3-2, as 
well as the absorption of the 12 CO 6-5 intensity. 



8 



Emprechtinger M. et al.: The Molecular Environment of NGC 2024 



The Fl component represents the radiation heated counter- 
part of B2 on the nearside of the HII region. This additional gas 
component does not change the modeled line emission as long 
as the H2 column density is lower than 4.6 • 10 2() cirT 2 , i.e. as 
long as its column density is smaller than 5 times the column 
density of B2. Because the existence of such a component is 
physically reasonable, we assumed that the interface between 
the bulk of the foreground material (F2) and the HII region is 
heated similar to the interface at the background, i.e., column 
density and maximum temperature of Fl and B2 are similar. 
Density and velocity of Fl are the same as found for F2. 

The geometric arrangement of the components Bl, B2, Fl, 
and F2 is well defined by our modeling and the temperature- 
and density trend of these components are shown in Fig. [5] 
However, the location of F3 is not clear, apart from residing 
in front of the HII region. The fact that this component is cold 
and not blueshifted with respect to the background component 
indicates that this component might be closest to the observer. 

A comparison between the observed spectra and the model 
results is shown in Fig. [4] Given the known complexity of the 
source and the limits of a five-component model, the model 
spectra fit the observations quite well. Two spectral features 
that were not included in the model fit are the redshifted out- 
flow, which can be seen in 12 CO 6-5, I2 CO 4-3 and 12 CO 3-2, 
and the component at Vlsr ~ 4.5 km/s. Nevertheless, the main 
features of the five emission lines can be explained with our 
rather simple assumptions. 

In Fig. [3] we show the integrated intensity of the 12 CO lines 
versus J up . Especially for the low- and mid-J CO lines the 
model matches the observations quite well. At J up > 10 the 
scatter of the observed data increases, possibly due to observa- 
tional uncertainties. 

5. Physical interpretation and discussion of the 
model results 

The model discussed above, successfully explaining the ob- 
served lines including the profiles of the velocity resolved lines, 
is the simplest model scenario satisfying the constraints im- 
posed by the blister and PDR picture. 

The Bl component is equivalent to the DMC in Fig. Q] 
Because the properties of Bl are mostly based on the fits of 
lines with J up < 5, only the warm gas of the DMC is repre- 
sented. Cold gas, which is deeper inside of the DMC, would 
have no effect on the modeled lines, either because the up- 
per levels are not populated ( 12 CO 13-12, 13 CO 6-5), or the 
lines are already optically thick ( 12 CO 7-6 and 12 CO 6-5). Thus, 
the column density of Bl is just the lower limit for the DMC. 
Observations of C ls O 2-1 and C n O 2-1 (Graf et al. [19931 re- 
vealed a three times lager H2 column density than what we 
found here. The B2 component represents the thin, hot interface 
between the HH-region and the molecular gas. Due to the high 
UV radiation emitted by the ionizing source(s) of NGC 2024, 
the gas at the interface is strongly heated, which may result in 
temperatures up to 330 K. 

The three subcomponents Fl, F2, and F3 together repre- 
sent the foreground material, t in Fig.Q] The Fl subcomponent 
is, as mentioned above, not traced by the lines we observed 
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Fig. 4. Observed (black) and model (red) spectra of 12 CO 13- 
12, 12 CO 6-5, 13 CO 6-5, 12 CO 3-2 and I3 CO 3-2. The observed 
spectrum of the 12 CO 13-12 line was shifted by -2.2 km/s (see 
text). 



here, but it represents the counterpart of B2 at the foreground 
component albeit at lower density. For the other two subcom- 
ponents located in front of the HII region, F2 and F3, we as- 
sume that matter is distributed uniformly across the area of the 
beam. Interstellar matter, however, shows clumpy structure on 
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Fig. 5. Temperature and density structure of NGC 2024. The 
observer is located to the right. The component F3 resides on 
the near side of the HII region, but its exact position is not 
clear and it might be physically disconnected from the other 
components. 



all scales. In NGC 2024, substructures on scales down to the 
resolution of previous observations were detected (Kramer et 
al.[l996 Lada et al. 119911 see SectionHJ. Therefore, we expect 
a clumpy structure within the beam as well, which would cause 
a much more complex structure of the foreground. Such small 
scale structures must be taken into account for more detailed 
modeling. For example, dense, cold clumps, smaller than the 
FWHM of the beam, which mainly affect the 12 CO 3-2 line, 
but due to their low temperature leave the I2 CO 6-5 line rather 
unaffected, are a possible explanation for the narrow absorption 
notches observed in 12 CO 3-2 and 4-3, whereas the line shape 
of 12 CO 6-5 is much smoother. This might be especially true 
for the thin and cold component F3. However, observations at 
higher spatial resolution would be necessary to verify this hy- 
pothesis. 

The velocity difference between the background gas and the 
bulk of the foreground material is 1.5 km/s. Assuming that the 
velocity difference is caused by the pressure of the ionized gas, 
according to the "Blister model" this velocity shift corresponds 
to a total momentum per unit area of 5670 g cirT 1 s _1 , which 
has to be transferred to the F2 component. We calculated the in- 
ternal pressure of the ionized gas, which pushes the foreground 
component, using measured electron temperature (8 160 K) and 
electron density of ~ 10 3 crrT 3 (Odegard 1985 ). Such a pres- 
sure has to act for approximately 2 • 10 5 years to accelerate 
the foreground component to its current velocity. However, the 
mass of the Fl component is highly unknown, and T e and n e 
may have changed in the past, so that this simple calculation 
may be somewhat debatable, but it shows that the observed ve- 
locity structure can be indeed explained by a blister model. 

5.1. Comparison with other models 

The properties of the background components agree fairly well 
with the results obtained by Graf et al. (I1993K The main dif- 
ference between their model and the present one is that Graf 
et al. dl993l l assumed LTE conditions, whereas in the SimLine 
code, the balance equations for all level populations and energy 



densities are solved self-consistently (Ossenkopf et al. |2001| l. 
In the previous study, an H2 column density of the background 
component of 2.0 ■ 10 23 cm 4 was found, which is three times 
more than what we found here. This difference arises because 
they included rarer isotopes, e.g., C ls O and C n O, which are 
optically thin and trace material located deep inside the cloud. 
Most of this material is hidden due to the high optical depths 
in the 12 CO and 13 CO lines we observed. The temperature de- 
termined by Graf et al. ( 11993I ), 67.4 K, is in good agreement 
with the 75 K determined in our work. Also the kinematic 
parameters Graf et al. (1993) report (vlsr = H-l km/s and 
Av = 1.8 km/s) match our results. Since they do not observe 
high rotational transitions of CO (no lines with J up > 7), they 
have no information about the hot interface region. 

Giannini et al. ( 120001 ) give a CO column density of 2 - 
5 • 10 18 cm 2 , which corresponds to an H2 column density of 
2.1 - 5.3 ■ 10 22 cirT 2 . This column density, which is about 
50% lower than the total column density we found (NH, ; totai — 
8.3 ■ 10 22 cirT 2 ), was derived using an escape probability code. 
The fact that Giannini et al. fitted only a single, isothermal 
gas component might explain the different results. In particular 
the background material will be underestimated in such an ap- 
proach due to self absorption effects. The kinetic temperature 
they give (110-130 K) is in the range of the values we found, 
but because of their isothermal approach hard to compare with 
our results. 

For the foreground components, we found a total H2 col- 
umn density of 1.07 • 10 22 cm 2 , which corresponds to an A v 
of 11.4. This is lower than the visual extinction of 15 to 25 
found towards several stellar objects inside the HII region (Bik 
et al. 2003). However, due to the low spatial resolution of our 
observations, the beam-averaged column density might be well 
lower than the values obtained from individual stars. The H2 
column density of the foreground component given by Graf et 
al. (119931) is a factor of two larger than our value. They also 
found a temperature of 23.5 K, i.e., slightly lower than our re- 
sults (30-40 K), which might be caused by the fact that they 
assumed LTE conditions. 

6. Conclusion 



We present the observations of seven 12 CO and 13 CO lines from 
J up = 3 to J up = 13 towards NGC 2024, a well know mas- 
sive star forming region in Orion. The shapes of these lines 
range from almost Gaussian ( 13 CO 6-5, 12 CO 13-12) to highly 
complex, multiple peaked ( 12 CO 3-2 , 12 CO 4-3), which indi- 
cates a complex internal structure of NGC 2024. In our anal- 
yses we also included the integrated intensities of 12 CO lines 
with J up > 14 (Giannini et al. 120001) . 

We modeled the high-J CO lines, which all seem to be 
optically thin, using an escape probability code (Stutzki & 
Winnewisser 1985) and found that these lines are emitted from 
a hot (> 250 K), dense (> 3 • 10 s cm 3 ), and thin (~ 100 AU) 
layer, which is most likely located at the interface of the HII 
region and the molecular cloud. 

We constructed a model for NGC 2024 using ID radiative 
transfer code SimLine (Ossenkopf et al. 2001). The velocity 
structure of this model is based on the principles of the Blister 
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model and the temperatures and column densities of the com- 
ponents are constrained by the PDR scenario. This model ex- 
plains the profiles of the observed 12 CO and 13 CO lines quite 
well. In our model, the bulk of the molecular gas resides at the 
back of the HII region and consists of warm (75 K) and dense 
(9 • 10 5 cm 3 ) material. We also find evidence for a hot (up to 
330 K) and thin (400 AU) layer located at the surface of the 
molecular cloud, from which most high-J CO emission origi- 
nates. The molecular ridge in front of the HII region consists of 
molecular material at lower densities (3 • 10 4 crrT 3 - 10 s cm -3 ). 
The column density is of the order of 10 22 cm 2 , which corre- 
spond to an A v of 1 1 mag. This value is in good agreement with 
measurements of the optical extinction towards stars within the 
HII region (Bik et al. 120051 . 

Overall, this study shows that for the example of 
NGC 2024, emission lines with complex and varying line 
shapes, as often observed in massive star forming regions, can 
be explained consistently with a rather simple, yet physically 
reasonable model. To explain the twelve emission lines a model 
of ~ 15 free parameter is required. Both the relatively large 
number of successfully fitted line intensities and line profiles 
and the fact that the multi-layer parameters are consistent with 
the physical scenario of a blister and a PDR, lead us to the con- 
clusion that we found a plausible model for the warm molecular 
gas in NGC 2024. 

Complex line profiles of low-J CO lines are commonly 
observed in massive star forming regions and with modern 
observatories (e.g., SOFIA, APEX and Herschel), more and 
more mid- and high-J CO observations will be available. This 
case study of NGC 2024 shows that with physical insight into 
these complex regions and careful modeling, the complex line- 
shapes and multi-line observations can be used to derive valu- 
able information on the physical conditions in massive star 
forming regions. 
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